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Abstract. 

We have carried out an extensive observational study (from BIMA data) and made a preliminary theoretical 
investigation of the molecular gas around HH 2. The molecular maps show a very complex morphological, kine- 
matical and chemical structure. For clarity we divided the observed region in four subregions: (1) The Ahead 
Core, located ahead of HH 2: its chemistry may be a consequence of a weak UV field originating in HH 2. The 
chemical structure within the Ahead Core suggests that it is not homogeneous but probably composed of small 
clumps; (2) The SO2 Clump, which is a molecular component within the Ahead Core that is more exposed to the 
UV radiation from HH 2. An increase of density and relative molecular abundances is observed towards HH 2. The 
UV radiation is possibly the source of molecular enhancement. Our chemical analysis confirms that this clump 
must have substructure within it; (3) the West Core, which is surrounded by a ring structure of shocked ionized 
gas and mid-IR emission. The ring structure is likely a consequence of the fact that the core is in the foreground 
with respect to the shocked and hot component. The chemistry of this core can be best explained as arising from 
a combination of an old photo-processed dense clump and a PDR, with or without a warm interface created in 
the interaction of the outflow with the core; (4) The High Velocity Region, associated with HH 2, is traced by 
HCO + but not by other molecular shock tracers. The chemistry can be accounted for by the interaction of the 
VLA 1 outflow with a dense clump via non-dissociative shocks and by the presence of a very strong UV field. The 
overall main conclusion of this work confirms the findings of Paper I and II, by demonstrating that in addition 
to the strong photochemical effects caused by penetration of the UV photons from HH 2 into molecular cloud, 
a range of complex radiative and dynamical interactions occur. Thus, despite the apparent 'quiescent' nature of 
the molecular cloud ahead of HH 2, the kinematical properties observed within the field of view suggest that it is 
possibly being driven out by the powerful winds from the VLA 1 protostar. 

Key words. ISM: individual: HH 2 — ISM: abundances — ISM: clouds — ISM: molecules — Radio lines: ISM 
— Stars: formation 



1. Introduction 

Very powerful winds are associated with the earliest stages 
of star formation. These winds alter significantly the 
dense molecular environments that surround the proto- 
stars. There are two kinds of interaction between the 
winds with the dense molecular environment; one is dy- 
namic, and the other radiative. The molecular outflows 
are the main signpost of the dynamical interaction (e.g., 
Richer et al. 12000(1 . The radiation generated in the strong 
shocks produced in the Herbig-Haro (HH) objects, as well 
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as the shock themselves alter the chemical composition 
of the molecular gas through mantle desorption and gas 
phase endothermic reactions (e. g., Girart et al. 1994, 2002 
Flower et al. 119961 Amin OTjOl]! . Thus, it is clear that high 
angular and spectral resolution maps of the molecular gas 
where this interaction takes place are very important in or- 
der to understand better the interaction mechanism and 
the properties of both the outflow and of the molecular 
clouds. 

HH 2 is a well studied bright, high excitation HH ob- 
ject, which generates strong UV radiation (e. g., Bohm et 
al. 119921 Raymond et al. I1997f) and presents a very com- 
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Table 1. Frequency setups of the BIMA observations 



Frequency (GHz) System 



LSB - 


- USB 


Main molecular transitions observed 




Dates 


Temperatures 


72.2 


- 75.8 


DCO+ 1-0, S0 2 6 ,6-5i,5, DCN 1-0, CH 3 SH 3 -2 A+ 




2000 Nov, 2001 Mar 


250-1000 K 


72.7 


- 76.4 


H 2 CO lo,i-0 ,o, S0 2 6o,6-5i, 5 , DCN 1-0 




2000 Apr 


300- 850 K 


85.3 


- 89.0 


HCN 1-0, HCO+ 1-0, H 13 CO+ 1-0, C 3 H 2 2i, 2 -lo,i, HCS' 


t" 2-1 


2000 Apr 


200- 600 K 


86.1 


- 90.0 


SO 2 2 -li, HCOOH 4 ,4-3 ,3 




2000 Jul 


200- 500 K 


86.7 


- 89.7 


H 13 CO+ 1-0, H 13 CN i-0, SiO 2-1 




2003 May, Jun a 


170- 600 K 


96.7 


-100.0 


CH 3 OH 2„-l„, HC 3 N 11-10 




1999 Oct 


200- 500 K 


98.0 


-101.0 


CS 2-1, OCS 8-7 




1999 Dec, 2000 May 


200- 800 K 


99.3 


-103.0 


SO 3 2 -2i, H 2 CS 3o,3 _ 2o,2 




2000 Jul 


300-1400 K 


100.5 


-104.1 


S0 2 3i, 3 -2o,2, NH 2 CN 5i, 4 -4i, 3 




2000 Nov, 2001 Mar, May 


140- 400 K 


104.2 


-107.8 


S0 2 3i,3-2o, 2 




2000 Jul 


300- 800 K 


109.7 


-113.0 


SO 2 3 -l 2 , C ls O 1-0, CN 1-0 




2000 Apr, May 


280-1300 K 


109.9 


-113.4 


C ls O 1-0, 13 CO 1-0, CN 1-0, CH 3 CN 6„-5„, HNCO 5 , 5 


-4 ,4 


2000 Nov, 2001 Mar-May 


350-1100 K 



a 2003 June observations were carried out in the D configuration 



plex morphology (Hester et al l 19981 Bally et al. l2002[l . It is 
driven by HH 1-2 VLA 1 (hereafter VLA 1), a Class pro- 
tostar (Andre, Ward-Thompson & Barsony|2000) . HH 2 is 
also associated with CO emission from the VLA 1 molec- 
ular outflow (Moro-Martin et al. 1993- A dense cloud 
of molecular gas and dust, apparently quiescent, appears 
ahead of HH 2 (Davis, Dent & Bell Burnell lTMll Torrelles 
et al.HnnS Dent, Furuya & Davis l2T)TI3| . 

Girart et al. (2002 hereafter Paper I) carried out a 
molecular line survey in the 0.8-4.0 mm range and at an 
angular resolution of ~ 30" of the dense molecular gas 
ahead of HH 2. From these observations we derive a tem- 
perature of ~ 13 K and a density of 3xl0 5 cm~ 3 . The 
observations show a characteristic chemistry that can be 
accounted for with chemical models of irradiated clumps 
(e.g., Taylor & Williams IT3§6l Viti & Williams HH)); the 
models indicate that photochemistry induced by HH ra- 
diation in regions that have not yet been shocked should 
produce a characteristic chemistry that is a signature of 
that origin. A detailed modeling of the conditions of the 
molecular clump ahead of HH 2 was carried out by Viti et 
al. H2003I hereafter Paper II). From the chemical analysis 
of this work we found that the 'illuminated' clump was 
young when it was first irradiated by HH 2 and that to 
account for the observed column densities a density gra- 
dient is required (i. e., the clump does not have a uniform 
density). In addition, we found that the peculiar chem- 
istry will last for no more than a few hundred years, since 
the UV radiation will ultimately destroy the molecular 
content. The chemical effects, however, last longer as the 
relevant region 'eats' into the core. 

In Paper I we studied the observational properties and 
the molecular gas at an angular resolution of 30" at one 
position ahead of HH 2, by using data from the CSO and 
BIMA. In Paper II we carried out chemical analysis by 
using the results of Paper I. The aim of this work is to 
study in detail, at higher angular resolution (~ 10"), the 



distribution and properties of the molecular gas around 
HH 2 and to understand how the powerful outflow arising 
from VLA 1 affects its properties. In § 2 we summarize the 
observational parameters of the maps. In § 3 we describe 
the morphological properties of the molecular emission. In 
§ 4 we analyze the data. In § 5 we discuss the chemical 
as well as the physical properties of the molecular envi- 
ronment around HH 2 and we consider why the regions 
observed appear to be chemically distinct. A brief sum- 
mary and the conclusions are given in § 6. 

2. Observations 

The observations were carried out with the 10-antenna 
BIMA array 1 at the Hat Creek Radio Observatory. The 
phase calibrators used were QSOs 0541—056, 0607—085 
and 0609—157. Absolute flux calibration was checked by 
observing Saturn or Mars. A total of 12 frequency se- 
tups were used, which covered frequencies between 72.0 
and 113.5 GHz. All the frequency setups were observed in 
the C configuration (i.e., antenna separation from 6.3 to 
100 m). At 3 mm this configuration provides a visibility 
coverage between 20 and 300 nanoseconds, i.e., an angu- 
lar resolution of ~ 9". Tabled shows for each frequency 
setup, the frequencies of the lower and upper side band 
(LSB and USB), the main molecular transitions observed, 
the observation dates and the range of system tempera- 
tures achieved. The digital correlator was configured to 
sample part of the 800 MHz wide IF passband in several 
windows, with adjustable frequency resolution. The typi- 
cal window for the line observations was configured with 
a 25 MHz bandwidth and 256 channels, giving a spectral 
resolution of 97.7 kHz. Thus, with the powerful correla- 
tor we were able to observe many molecular transitions. 

1 The BIMA array is operated by the Berkeley-Illinois- 
Maryland Association with support from the National Science 
Foundation. 
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Table 2. Parameters of the BIMA maps 



Synthesized 







Beam 




rms 


Molecule 


V 


HPBW, PA 


Av 


(Jy/ 


Transition 


(GHz) 


(arcsec, deg) 


(km s ) 


beam) 


DCO + 1-0 


72.0393 


19.3x7.3, 1 


0.20 


0.30 


JtloC^CJ In i — On n 


72.8380 


18.0x9.4, 15 


0.20 


0.27 


C^^ri9 zi 9 — In i 


85.3389 


14.1x7.1, -1 


0.34 


0.10 


SO 2 2 -li 


86.0934 


12.9x7.6, 


0.34 


0.16 


H 13 CN 1-0 


86.3402 


17.6x11.8, 4 


0.33 


0.11 


H 13 CO+ 1-0 


86.7543 


17.4x11.2, 9 


0.33 


0.11 


SiO 2-1 


86.8470 


17.7x11.7, 6 


0.33 


0.11 


HON 1-0 


88.6318 a 


12.7x7.1, 


0.33 


0.10 


HCO+ 1-0 


89.1885 


13.2x7.0, 1 


0.33 


0.11 


CH3OH 2 n — l n 


96.7414 6 


11.7x6.4, 4 


0.30 


0.13 


CS 2-1 


97.9810 


14.2x6.1, 10 


0.30 


0.14 


SO 3 2 -2i 


99.2999 


11.9x6.6, 15 


0.59 


0.14 


SO2 3l,3-2o,2 


104.0294 


9.9x5.6, 1 


0.28 


0.05 


SO 2 3 -l 2 


109.2522 


10.4x5.9, 2 


0.27 


0.15 


c 18 o 1-0 


109.7822 


10.1x5.9, 1 


0.27 


0.12 


13 co 1-0 


110.2014 


9.7x6.3, -8 


0.27 


0.13 


CN 1-0 


113.4910 c 


12.0x8.9, -5 


0.52 


0.08 



a Frequency for the F=2-l hyperfine line. Observations also 
detected the F=l-1 and F=0-1 hyperfine lines 

6 Frequency for the 2q-1oA + line. Observations also detected 
the 2_i-l_i E line 

c Frequency for the 1-0 §, §-§, § line. Observations also de- 
tected the §, 5—51 \ and |, 5—5, § lines. The maps were 
obtained by applying a Gaussian taper of 26 kA to the vis- 
ibilities. 



Table 2 from Paper I gives the list of the lines observed 
before 2003, with the peak intensity or upper limits for an 
angular resolution of 30". In the 2003 observations, of the 
three lines observed (see Table we only detected the 
H 13 CO+ 1-0. 

The phase center of the observations was a( J2000) = 
5 h 36 m 27?20; 5(J2000) = -6°47'27". The FWHM BIMA 
primary beam ranges from 2'.7 at 72 GHz to 1'.7 at 
113 GHz, and thus our observations engulfs within the 
primary beam the HH 2 object and dense core associated 
with it. For the 109.9-113.4 GHz (LSB-USB) frequency 
setup, the observations (which included the 13 CO, C 18 
and the CN lines) were done with a mosaic of five points: 
a field centered in the phase center, and the other fields 
located 48" N, S, E and W of the the phase center. The 
data were calibrated using the MIRIAD software. Maps 
were made with the (it, v) data weighted by the associated 
system temperatures and using natural weighting. Tabled 
lists all the molecular lines detected from the observations 
listed in Paper I as well as the three lines observed in 
2003. The table includes the channel resolution, the re- 
sulting synthesized beam and the rms noise of the maps 
at this channel resolution. 




Fig. 4. Superposition of the [SII] image (from Curiel, pri- 
vate communication) and the BIMA contour maps of the 
zero-order moment (integrated emission over the 5.3 to 
11.3 km s" 1 velocity range) of the HCO+ 1-0 and 13 CO 
1-0 around HH 2. The BIMA maps are corrected by the 
primary beam response. Contours are 1, 2.5, 4, . . . , 32.5 
times 0.45 (HCO+) and 0.72 ( 13 CO) Jy km s" 1 . The con- 
tours with dots around HH 2 indicates an emission valley. 
The dashed thick line shows the BIMA beam response at 
the 0.5 level with respect to phase center. The synthesized 
beam is shown in the bottom left corner of the panels. 
The position of the driving source of the HH 1-2 outflow 
(VLA 1) is shown as a filled square in the upper-right 
corner of the panels. 

3. Results 

Figures ^ HI an d B1 show the channel maps over the 6 to 
10 km s _1 velocity range of the molecular transitions 
with emission detected at least in one channel. For most of 
the molecules detected their emission is comprised within 
this velocity range. Figure^shows the zero-order moment 
of the emission for the 13 CO and HCO+ 1-0 lines. The 
kinematical and morphological complexity of this region 
is clear from these figures. 

For clarity, we will divide the observed region in four 
subregions taking into account the properties of the dif- 
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Fig. 5. Sketch of the HH 2 region. Thick grey contours 
show the main molecular regions (the Ahead and West 
Cores with respect to HH 2 and the VLA 1 dense core). 
The dashed contours show the high velocity HCO + clump. 
The dotted points contours show the jet-like molecular 
outflow (obtained from the HCO+ emission at a wlsr ve- 
locity of 7.3 km s _1 after masking the contribution from 
the Ahead Core). The grey scale shows the [SII] emission 
(from Curiel, private communication). The white labels 
within the [SII] emission indicates the HH 2 knot names. 
The small filled ellipsoid shows the position and elonga- 
tion direction of the VLA 1 radio jet, the powering source 
of the HH 1-2 outflow. The arrow shows the direction of 
the outflow. 



ferent molecular species detected and the properties of 
other tracers detected at other wavelengths (see Figure|5J): 
(1) The molecular core ahead of HH 2 (hereafter Ahead 
Core); this region includes the region studied in Paper I 
and II, which coincides with the compact emission traced 
by SO and SO2; (2) The molecular core west of HH 2 
(hereafter West Core); (3) The high velocity molecular 
emission spatially associated with HH 2 (hereafter High 
Velocity Region) ; (4) The dense core surrounding the HH 2 
driving source, VLA 1. In the following three subsections 
we describe the first three regions. The VLA 1 dense core 
is detected clearly in several molecules ( 13 CO, H 13 CO + , 
H2CO, HCN, C3H2) in the 10 km s _1 channel map (see 
Figs, n [21 an d EJ- However, since it is located far from 
the phase center of the BIMA observations (~ 1'.7) we do 
not study its properties. The properties of this core have 
been well described in previous studies (e.g., Torrefies et 
al.HnHI Choi & Zhou [TM7| . Thus, below we describe the 
first three regions. We also describe the overall kinematics 
of the dense gas seen within the field of view. In the last 



subsection of the results we briefly describe the emission 
traced by the 13 CO 1-0. 

3.1. The Ahead Core 

We define the Ahead Core as the cold dense molecular and 
dust component located ahead of HH 2; it was first de- 
tected from observations of HCO + (Davis et al. H990J> and 
NH 3 (Torrelles et al.\M2\MB- Dent et al. (|2UU5|l have 
recently studied the submm dust emission from the Ahead 
Core, which arises from a region with a size of 0.13x0.09 
pc 2 . From the dust properties, Dent et al. ( 2003J) derive 
a total mass for the core of 3.8 Mq, an averaged dust 
temperature of 22 K and a core luminosity of 13 Lq. 

For some species (e.g., HCO + : see Fig. 0J, the maps 
of their integrated emission over all the velocity range 
where the emission is detected (excluding the high ve- 
locity component) show that the Ahead Core is not an 
isolated molecular structure but it is connected with the 
VLA 1 dense core. Indeed and despite the BIMA primary 
beam attenuation, the overall morphology of the HCO + 
1-0 agrees well with that from the HCO + 3-2 integrated 
emission (Choi & Zhou C 997). However, we consider the 
Ahead Core as the molecular emission that is spatially co- 
incident with the submm dust structure detected ahead of 
HH 2 (see Figure 2 of Dent et al. 2003). The emission as- 
sociated with the Ahead Core emits in the 5.2-7.8 km s _1 
wlsr range. Figure shows the averaged emission within 
this velocity range for some of the molecular transitions 
presented in this paper as well as NH3 (1,1) from Torrelles 
et al. I|1994fl . and the 0.85 mm dust map from Dent et al. 

mm . 

The emission of the Ahead Core is traced by all the de- 
tected lines in this work, except for the C3H2 2^2-10,1 line 
(which is only detected towards the VLA 1 dense core), 
as can be seen in the 6, 7 and 8 km s _1 channel maps of 
Figures n El an d El Despite the apparent complexity of the 
molecular emission for the different species, we can group 
them according to their morphologies within the Ahead 
Core. 

First, HCO + and isotopes, as well as HCN and H2CO 
are the only species that trace the whole dusty structure 
from the Ahead Core. In particular, Fig. H3 shows that 
there is a good spatial correspondence of these species 
with the submm dust emission. NH3 can also be included 
in this group, since its emission also follows most of the 
dust emission (Fig. EJ), especially if compared with the 
molecular emission from species of other groups (see be- 
low). In spite of the overall agreement between the emis- 
sion of these molecules and the dust, the intensity peak 
of the molecular emission is located around HH 2 knot L 
(only HCO + peak coincides with that of the dust). HCO + 
is the only molecule that has emission engulfing most of 
HH 2, including the brightest and highest excitation knots, 
H and A, where no dust is detected. 

A second group of molecules comprise SO, SO2 and 
CH3OH: they appear to be significantly more compact 




Fig. 6. Contour maps of the ahead core for the dust emission at 0.85 mm (adapted from Figure 2 of Dent et al. 
2003) and for several molecular lines. The grey scale [SII] image (from Curiel, private communication) is also shown 
in the molecular line panels. The molecular maps shown here are obtained by averaging the emission over the 5.2 to 
7.8 km s _1 wlsr velocity range. The NH3 (1,1) map was obtained from the data of Torrefies et al. ( 1994). The triangle 
shows the position of the 0.85 mm dust emission peak. The dotted rectangle in the HCN panel shows the area used 
to compute the molecular abundances in the Ahead Core. 



than the previous species (and the dust) and their emis- 
sion is located ahead of the brightest HH 2 knots and with 
the strongest emission spatially coincident with HH 2 knot 
L (Fig. EJ. Interestingly, the emission of CH3OH, SO and 
SO2 appears displaced to the west with respect to dust 
peak, i. e., apparently, they trace the face of the dust core 
exposed to the HH 2 object (see Fig. |BJ). 

A third group includes CS and the CO isotopes. The 
CS and C 18 emission and the bulk of the 13 CO emission 
from the Ahead Core arise south of the HH 2 knot L, with 
little or no emission in the east part of the Ahead Core. 
The peaks of emission of these molecules are located 15" 
to 30" south of the dust peak and are roughly coincident 
with secondary dust peaks. 

Finally, the CN emission has an inverted 'L' morphol- 
ogy centered and peaking at the same position as the dust. 
To the south of the dust peak it follows approximately the 
CS emission. 



(Dent et al. I2003|> . Interestingly, very recent mid-IR ob- 
servations shows that this molecular clump is surrounded 
by a ring structure of hot dust (Lefloch et al. 2005). This 
hot ring shows also PAH emission (Lefloch et al. 2005), 
Ha (Warren-Smith & Scarrott H^jl and [SII] (Fig.Q. 

The West Core is detected in HCO+, 13 CO, C 18 0, 
H 2 CO, CN, HCN, CS and marginally also in S0 2 and 
C3H2. This core is not detected in NH3 (Torrefies et al. 
1994). Figure [7| shows clearly how well the [SII] surrounds 
the molecular emission. The [SII] hole, where the molecu- 
lar emission appears, is remarkably circular with a diam- 
eter of ~ 30" or 0.06 pc. The [SII] hole is especially well 
traced by the CO isotopes, whereas the CN and HCO + 
are slightly shifted to the north with respect to the center 
of this hole. The molecular line widths are quite narrow 
(Av ~ 0.4 km s _1 ) and centered at ulsr of 6.3 km s~ x . 
A weaker and broader component is also detected at vlsr 
= 8.8 km s _1 . 



3.2. The West Core 

The West Core is the molecular structure detected west 
of HH 2 at the 6 km s" 1 channel of Figs. [T] H and [2 
This core is also traced by the dust at submm wavelengths 



3.3. High Velocity Region 

Similarly to Dent et al. QfflEft, we found broad HCO+ 1-0 
emission (i>lsr ranging from —8 to 24 km s^ 1 ) spatially 
compact and associated with HH 2. Figure|S]shows clearly 
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Fig. 7. Bottom panels: Superposition of the grey scale 
[SII] image (from Curiel, private communication) and the 
13 CO, HCO+, H 2 CO, C 18 0, CN and HCN BIMA con- 
tour maps at the 6 km s _1 velocity channel of the 
molecular clump west of HH 2. The [SII] image was con- 
volved with a Gaussian with a FWHM of 4" in order 
to enhance the weak ring structure around the molecu- 
lar emission. Contours are 50, 70 and 90% of the peak. 
The cross shows the position of the submm dust posi- 
tion shown in Figure 2 of Dent et al. i2003[). The syn- 
thesized beam of the molecular maps are shown in the 
bottom right corner. The dashed circle shows the area 
used to estimate the column densities (see § I4.2f> . Top 
panels: Primary beam corrected spectra of the molecu- 
lar lines (shown in the bottom right of the panels) taken 
at a(J2000) = 5 h 36 m 23?2; <5(J2000) = -6°47'30" and at 
an angular resolution of 30". For the CN the transition 
shown is the N=l-0 J=3/2-l/2 F=5/2-3/2. The 13 CO 
and HCO + line intensities have been scaled down a factor 
3 and 2, respectively. 



that the high velocity HCO + emission is well correlated 
with the shock-excited near-IR l-O S(l)H2 emission. The 
spatial correlation is even better with respect to the H 2 
0-0 S(2) line (see Fig. I from Lefloch et al. ED03J: The 
strongest high velocity HCO + emission appears located 
around the HH 2 knot E and there is also weaker emission 
near the highest excitation HH 2 knots. A more detailed 
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Fig. 8. Bottom panels: Superposition of the grey scale 
near-IRH 2 2 1-0 S(l) image (K' band; Curiel, private com- 
munication) and the contour map of the integrated high 
velocity HCO + emission over the —5 to 5 and 11 to 
21 km s" 1 i>lsr interval, respectively. HCO + contours 
levels are 3, 5, 8, 12, 17, 22, 27, 32, 37 and 42 times the 
rms noise of the map, 0.2 Jy km s" 1 . The synthesized 
beam is shown in the bottom left corner of each panel. 
Top panel: HCO + 1-0 spectra of the high velocity gas. 
The cross marks the position of the X-ray source (Pravdo 

et ai.nnnj. 
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Fig. 9. Spectra of the different molecular transitions at 
the intensity peak of the high velocity HCO + 1-0 emission 
(a(J2000) = 5 h 36 m 26!57; <5(J2000) = -6°47'9'.'l). 
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description of the high velocity HCO + around HH 2 and 
its properties is given by Dent et al. ( 2003 ) . 

As shown in Figure |5| and from our BIMA data, 
the high velocity emission associated with HH 2 is de- 
tected only in the HCO + 1-0 line (the only other re- 
ported molecule with high- velocity emission is CO: e.g., 
Moro-Martin et al. 11999(1 . This is interesting, since species 
that are usually strongly enhanced in shocked molecular 
clumps, such as SiO, CS, CH 3 OH , H 2 C O and SO (e.g., 
L1157: Bachiller & Perez Gutierrez ll997f) are not detected 
in this region. 

3.4. The overall kinematics of the dense gas 

A look at the channel maps of the HCO + emission (Fig.[IJ 
shows that there is a clear velocity gradient in the dense 
molecular structure connecting the Ahead Core with the 
VLA 1 dense core. That is, the emission closer to the sys- 
temic velocity (wlsr — 9.5 km s , Moro-Martin et al. 
1999) appears located closer to the VLA 1 dense core, 
whereas the emission with bluer velocities appears located 
closer to the Ahead Core. This is also observed in the emis- 
sion of other molecules but with a weaker signal-to-noise 
ratio (e. g., HCN and H 2 CO). Note that the emission from 
this molecular structure from the 9 and 10 km s _1 velocity 
channels of Figures ^ 121 an d CD is different for the different 
molecules. Since this is the systemic velocity, these dif- 
ferences could be due to high optical depths, or to other 
factors, such a different excitation or chemical conditions. 
We will not discuss further these differences and we will 
discuss only the kinematics from the HCO + and H 13 CO + . 

The velocity gradient is also clearly shown in the first- 
order moment map of the HCO + (Fig. 110(1 . The most 
blueshifted velocities, wlsr ~ 6.4 km s _1 , appear close 
to HH 2 knot L in the Ahead Core, and in the West Core. 
South of the Ahead Core the velocity of the gas goes back 
to the systemic velocity. 

The second-order moment map of the HCO + (Fig. 110(1 
shows that most of the emission has a line width (FWHM) 
in the 0.7-1.6 km s _1 range. The largest line widths are 
observed around HH 2 knot E, where the high velocity 
HCO + is observed. Large line widths are also observed 
along three arms emerging around knot E: one connecting 
with VLA 1 and coinciding with the outflow axis and the 
other two perpendicular to this one. This region of line 
width enhancement form a rotated 'T-like' shape. 

In order to see better the overall kinematics and, in 
particular, the region of line width enhancement, we show 
in Figure El the position-velocity maps for the HCO + 
and H 13 CO + 1-0 along three different cuts. Two of them 
(cuts 1 and 2) are along the direction of the line widths 
enhancement (see previous paragraph). The other (cut 3) 
is parallel to cut 2 but further away from VLA 1 (although 
still within the Ahead Core). 

Despite of the complexity of Figures an d El the 
kinematics of the region can be summarized as follows: 




Fig. 10. Superposition of the the BIMA contour maps 
of the zero-order moment (integrated emission) with 
the false-color image of first-order (bottom panel) and 
second-order (top panel) maps of the HCO + 1-0. The 
moment maps were obtained in the 5.35-11.26 km s 
«lsr range. The cross shows the position of the 0.85 mm 
dust emission peak. The labels indicate the HH 2 knot 
names. 



— The region around HH 2 knot E shows the largest line 
widths because of the contribution from the high ve- 
locity HCO + (this corresponds to the 0" offset position 
in the 1 and 2 cuts of Fig [TT|l . 

— One of the regions of line width enhancement (the 
one connecting VLA 1 with HH 2) is spatially co- 
incident with the molecular outflow associated with 
HH 2. The HCO + emission presents a double peak 
along the molecular outflow axis (negative offset po- 
sitions in cut 1 of Fig. Ill|) . The emission from the 
more blueshifted peak arise from a collimated struc- 
ture connecting VLA 1 with the Ahead Core (see Fig. El 
and the 7 km s" 1 channel of the HCO+ of Fig. [TJ. 
Interestingly, the CO emission from the molecular out- 
flow in this region is mainly redshifted with respect to 
the systemic velocity (see Moro-Martin et al. [1999), 
which is the reverse situation of the HCO + emission. 
Taking into account that the outflow axis is almost on 
the plane of the sky, this suggests that the HCO + is 
possibly tracing the surface of interaction of the molec- 
ular outflow with the dense molecular gas in the fore- 
ground face of the outflow lobe. 

— Cut 3 shows that the most blueshifted velocities ap- 
pear at the intersection of the outflow axis with the 
Ahead Core (0" offset position) and the velocity be- 
comes closer to the systemic velocity as the emission 
gets further away from the outflow axis. The same be- 
havior is seen in the velocity gradient shown by the 
first-order moment map of Figure Cut 2 shows a 
similar trend than cut 3 but, in addition, there are two 
other velocity components that appears only in HCO + , 
the high velocity component at the 0" offset position 
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<1 




tot (") 

Fig. 11. Bottom panel: Zero-order moment map of the 
HCO + . The dashed lines and numbers shows the direction 
of the position-velocity cuts shown in the upper panels. 
The cross shows the position of the 0.85 mm dust emission 
peak. The labels indicates the HH 2 knot names. Upper 
panels: Position-velocity maps of the HCO + and H 13 CO + 
1-0 emission. The origin of position offsets of panels 1 
and 2 are the intersection of the dashed lines 1 and 2 in 
the bottom panel, while for panel 3 is the intersection of 
dashed lines 1 and 3. The sign of the offsets corresponds 
to the sign of Aa in the bottom panel. 

and a component at ~ 8 km s . Since cut 2 follows 
the border of the Ahead Core facing VLA 1 , the latter 
component may be an indication of interaction of the 
outflow with the dense core. 



3.5. The 13 CO emission 

The 13 CO emission (0 and Fig. shows significantly 
more structure than the other molecular species, includ- 



ing HCO + . This is because of the CO low dipole moment, 
which allows the molecule to trace low density molecular 
gas (undetected with the higher density tracers). Fig. 0] 
shows that within the 13 CO field of view (which is larger 
than for the other molecules because of the mosaiced ob- 
servations) the emission spreads out almost everywhere 
except in a small region spatially coincident with HH 2. 
The strongest emission surrounds HH 2. It is interesting 
to note that some of the structures traced by the high 
density tracers are not well traced or traced significantly 
differently by the 13 CO. 

At the systemic velocities (9 and 10 km s _1 velocity 
channel maps of Fig.^J) the strongest 13 CO emission arises 
from a clear V-shaped structure, with its vertex located 
near VLA 1, facing and surrounding HH 2. The northern 
part of this structure corresponds approximately to the 
high density molecular gas that connects the Ahead Core 
with the VLA 1 dense core. Given the direction of the 
outflow (from VLA 1 to HH 2), the V-shaped morphology 
seems to trace the walls of a cavity created by the HH 1-2 
outflow. Indeed, within this cavity and connecting VLA 1 
and HH 2 there is a highly collimatcd molecular outflow 
(Moro-Marti'n et al. H555^ . 

4. Analysis 

In order to compare properly the different species and de- 
rive the temperature and column densities, we first made 
maps of the different transitions at the same angular 
resolution. This angular resolution was set to 15" x 9", 
PA ~ 0°. 



4.1. The Ahead Core 

The physical parameters for the Ahead Core were derived 
from the maps of the integrated emission within the 5.2— 
7.8 km s _1 velocity range. 

4.1.1. Optical depth of the HCO+ and HCN 1-0 lines 

The HCO + 1-0 optical depth was derived from the ratio 
of the H 13 CO + and HCO + integrated emission maps, cor- 
rected for the primary beam response. The highest values 
of the HCO+ 1-0 optical depth within the FWHM of the 
primary beam are found south of HH 2, where the opti- 
cal depth reaches up to ~ 20. The optical depth smoothly 
decreases to 2-4 east and northeast of HH 2. 

The H 13 CN 1-0 is not detected, so only upper lim- 
its can be derived for the HCN 1-0 optical depth. Since 
the HCN is quite extended, channel maps of both isotopes 
were smoothed to an angular resolution of 20". From these 
maps we estimate that thcn £5 10. This implies that the 
observed HCN column density derived in Paper I could be 
underestimated by a factor of SS 10. In fact, from Paper 
II, the best matching model (B15) at ~ 100 yrs gave HCN 
abundances ranging from ~ 2xlO n cm~ 2 at 1 mag to ~ 
6xl0 16 cm -2 at 5 mags; for a 3 mags gas we have a the- 
oretical column density of 2.7xl0 13 which is indeed only 
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Table 3. Mean column densities and abundances in the 
Ahead Core 



-B*47' - 



Molecule 


JV[mol] 
( cm" 2 ) 


X[mol 


6 


ccr 


2.9xl0 17 


3.2x10" 


-5 


H 2 CO 


3.0xl0 13 


3.4x10" 


-9 


HCO+ 


2.9xl0 13 


3.2x10" 


-9 


CH 3 OH 


1.4xl0 13 


1.5x10" 


-9 


SO 


6.7xl0 12 


7.4x10" 


-10 


HON 


1.5xl0 12 


1.6x10" 


-10 


SO2 


1.2xl0 12 


1.3x10" 


-10 


CS 


9.2x10" 


1.0x10" 


-10 


CN 


6.7x10" 


7.4x10" 


-11 


DCO+ 


2.2x10" 


2.4x10" 


-11 



Here after the CO column densities are derived adopting a 
12 CO to 13 CO ratio of 63 (Langer & Penzias 115531 . 
X[mol] is the abundance with respect to H2 



a factor of 10 higher than the observed column density of 
4xl0 12 cm- 2 . 

4.1.2. Mean physical parameters 

The total mass of the Ahead Core is 3.8±0.4 Mq, as de- 
rived from submm dust observations (Dent et al. l2003|l . To 
estimate the mean column densities and molecular abun- 
dances of the observed species over the whole Ahead Core 
(using the mass derived from the dust), we integrated the 
emission of the different species within the 5.2-7.8 km s _1 
velocity range and over a box of 1 '0 x 1.'5 that encom- 
passes the dust emission ahead of HH 2 shown by Dent et 
al. I|2003|l . The column densities were estimated assuming 
an excitation temperature of 10 K. Table [3] shows the es- 
timated absolute molecular abundances averaged over the 
whole Ahead Core. The relative abundances with respect 
to the CO are a factor ~2-10 lower than the derived from 
Paper I. Note however that the relative abundances given 
in Paper I were obtained for a FWHM beam of 30" at 
the SO peak, whereas the values derived in Table were 
obtained over a larger region. This difference can be in- 
terpreted as arising from chemical gradients within the 
Ahead Core. 

4.1.3. Gradients within the Ahead Core 

We computed the excitation temperature of the SO for 
the 32~2i and 23-I2 lines. By using these two transitions 
the excitation temperature can be expressed as: 



11.7 log(1.77 



/ J(SQ2 3 -l 2 )<fo. 
//(SO 3 a -2t)dv 



(1) 



where / is the line intensity (in K). The resulting excita- 
tion temperature map is shown in Figure 1121 there is a 



JO" 



30" 




28' 



Fig. 12. Top panel: Superposition of the 

10- 5 A[SO 2 ]/A[CO] map (thick grey contours), the 
zero-order moment map of the SO2 3i,3-2o,2 emission 
(dashed contours) and [SII] emission (grey scale image). 
Bottom panel: Superposition of the temperature map 
(thick grey contours) derived from the SO 2 3 -l 2 to SO 
3 2 -2i ratio, the zero-order moment map of the SO 3 2 -2i 
emission (dashed contours) and the [SII] emission (grey 
scale image). The dotted line shows the slice done for 
Figure ED 



gradient, with increasing excitation temperatures towards 
HH 2. 

Column density maps for SO2 and 13 CO were derived 
assuming an excitation temperature of 10 K. Figure 1121 
shows the map of the SO2 relative abundance with re- 
spect to the CO. The map shows a gradient in the S0 2 
relative abundance, increasing towards HH 2, having the 
highest value close to the HH 2 knot E. Note that, if in- 
stead of using a constant T cx , we use the values derived 
from the SO line ratio analysis (equation 1) the results do 
not significantly change. 

Given the SO excitation temperature and SO2 relative 
abundance gradient towards HH 2 shown in Figure^] we 
measured the integrated emission along the direction of 
this gradient, which approximately coincides with the SO 
and SO2 major axis. From these slices we computed: (1) 
The excitation temperature for the SO 32~2i and 23-I2 
lines (see first paragraph of this section). (2) The column 
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Fig. 13. Slices of the SO excitation temperature (top 
panel), column density (central panel) and relative abun- 
dance with respect to the CO (bottom panel) as a function 
of the distance to HH 2. The first position in the plot coin- 
cides with the position of the HH 2 knot E. The slice has 
position a(J2000) = 5 h 36 m 26!5; S( J2000) = -6°47'06".'0 
and the position angle of the slice is PA=169°. Note that 
the difference in column density with Figure 1 from Paper 
II are because in paper II we used a constant T ex whereas 
in this paper we used the value derived from the SO line ra- 
tio. In the case of HCO + , the new observations of H 13 CO + 
allowed to correct for optical depth. 



density of the different molecular species assuming an ex- 
citation temperature equal to the SO excitation tempera- 
ture. For HCO + , its column density was corrected by tak- 
ing into account thco+ wherever H 13 CO + was detected, 
otherwise t HC q+ = 6.0 was adopted (which is the low- 
est Thco+ value measured in the slice). (3) The relative 
abundances with respect to CO. 

Figure 1131 shows the three parameters, SO excitation 
temperature, column density and abundance, along the 
slice. Table 0] gives the SO excitation temperature, col- 
umn density at four positions along the slice, as well as the 
abundance increase factor at the first and last positions. 



Table 4. Column densities for selected positions along the 
slice in the SO2 Clump 



+" 1 An 


h- 1 


12 


13 


14 




Tex (W 


27±9 


16±2 


13±2 


9±3 




rfslice (") b 


10 


20 


30 


40 




d HH2 cr 


10 


22 


29 


39 








N 


mol] 






Molecule 




( cm" 2 ) 






CO 


1.7xl0 17 


3.2xl0 17 


3.3xl0 17 


3.2xl0 17 


1 


H2CO 


l.OxlO 14 


1.2xl0 14 


9.0xl0 13 


4.8xl0 13 


4±1 


HCO+ 


6.7xl0 13 


6.0xl0 13 


4.0xl0 13 


1.8xl0 13 


8±1 


CH3OH 


8.0xl0 13 


7.1 xlO 13 


6.2xl0 13 


2.8xl0 13 


6±2 


SO 


1.6xl0 13 


2.1 xlO 13 


1.7xl0 13 


5.7xl0 12 


5±1 


HCN 


3.9xl0 12 


3.7xl0 12 


2.8xl0 12 


1.5xl0 12 


5±1 


S0 2 


1.4xl0 13 


9.7xl0 12 


6.0xl0 12 


<2xl0 12 


£12 


CS 


9.1x10" 


1.8xl0 12 


2.0xl0 12 


2.0xl0 12 


0.8±0.3 


CN 


4.1xl0 12 


3.7xl0 12 


3.0xl0 12 


2.4xl0 12 


3±2 


DCO+ 


1.5xl0 12 


9.9x10" 


7.8x10" 


<5xl0 11 


£ 5 



a Error bars are at 1-a 

b Relative distance in the slice (see Fig. I Kit 
c Distance to HH 2 knot H 



It is clear from this figure that the relative abundance 
enhancement towards HH 2 is observed in all molecular 
species except CS, which shows a relative abundance al- 
most constant along the slice. HCO + and SO2 are the 
species that have the highest enhancement. 

In order to interpret the excitation temperature en- 
hancement we used the RADEX package to fit the ob- 
served SO 32-2 1 / 2 3 -l 2 line ratios at the positions II 
(2.7±0.4) and 13 (4.4±0.6) given in Table H RADEX is 
a non-LTE molecular radiative transfer in an isothermal 
homogeneous medium (Schoier et al. I2005[) . In Figure PHI 
we show a set of the RADEX solutions in the Tki n -n(H2) 
plane for the observed line ratios. The set of solutions in 
Fig. [21 that satisfy these averaged values indicates that 
the density in II (il 10 6 cm -3 ) is higher than in 13 
(~ 10 5 cm -3 ), although the temperature cannot be con- 
strained. Further observations of higher excitation lines at 
high angular resolution will help to better constrain the 
density enhancement and to determine the temperature. 

4.2. West Core 

We first corrected the channel maps from the BIMA pri- 
mary beam response and then convolved them with a 
Gaussian with a FWHM of 30" in order to optimize 
the signal-to-noise ratio. The column densities were esti- 
mated from the resulting map at the position of the dust 
peak (see Figure[7J) and for the i>lsr of 6.3 km s _1 compo- 
nent. Since we do not have information about the temper- 
ature of the molecular gas (no SO line is detected towards 
this core), we assumed an excitation temperature of 10 K. 
The narrow line widths of this velocity component indi- 
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Table 5. 

Core 



Column densities and abundances of the West 





AT[mol] 
( cm- 2 ) 


X[mol]/ 






Molecule 


X[CO] 


X[mol] 


a 


CO 


1.7x10" 


1.0 


1.8x10 


— £> 


H 2 CO 


1.3xl0 13 


8xl0" 5 


1.4x10 


-9 


HCO+ 


6.8xl0 12 


4xl0" 5 


7.2x10" 


10 


CH 3 OH 


^2xl0 13 


;$lxlO" 4 


^2x10 


-9 


SO 


^2xl0 12 


;$lxlO" 5 


;$2xl0" 


10 


HCN 


1.6xl0 12 


9xl0" 6 


1.7x10" 


10 


S0 2 


~2xl0 12 


lxlO" 5 


~2xl0" 


10 


CS 


1.6xl0 12 


9xl0" 6 


1.7x10" 


10 


CN 


3.3xl0 12 


2xl0" 5 


3.5x10" 


10 


DCO+ 


~2xl0" 


lxlO" 6 


~2xl0" 


11 


C3H2 


~7xlO n 


4xl0" 6 


~7xl0" 


11 



X[mol] is the abundance with respect to H2 



log[n(H 2 )/cm 3 ] 

Fig. 14. Set of RADEX solutions in the Tki n -n(H2) plane, 
for the observed SO 32-2i / 23-I2 line intensity ratios 
(assuming a SO column density of 2xl0 14 cm" 2 ). The 
thick dashed lines is the contour for the ratio at position 
13 (4.4) and the thick dotted-dashed line is the contour 
for the ratio at position II (2.7). The grey area shows the 
area within a ±ler. 



cate that the molecular gas is not shocked and, therefore, 
the temperature is probably low: even in the case that 
the West Core is in a similar situation as the SO2 Clump 
(see § I5.3fl . the region with high temperatures would be 
probably compact. Fig. 2 of Dent et al. (20031) shows that 
the 850 /im dust emission peak towards the West Core is 
~ 90 mJy beam -1 . Assuming that the dust temperature 
and opacity is similar to the Ahead Core, then the mean 
gas column density in the West Core at the dust peak is 
iV(H 2 ) ~ 1 x 10 22 cm -2 . The beam averaged 13 CO column 
density at the dust peak and at an 14" angular resolution 
(the same as the dust map) is 2.8xl0 15 cm -2 , which 
gives an abundance of A[CO]~ 1.8xl0" 5 . Table shows 
the derived column densities and relative abundances for 
the 30" smoothed spectra (adopting the aforementioned 
X[CO] value). From the 13 CO emission, we derive a mass 
of - 0.5 M Q for the West Core. 



Table 6. Column densities and abundances in the High 
Velocity Region 



NCmolO/NCHCO" 1 





jV(mol) 




Shocked 


Model 


Model 


Molecule 


( cm" 2 ) 


HH 2 


Regions" 


A b 


B" 


HCO+ 


4.1xl0 13 


1 


1 


1 


1 


CO 


<lxl0 17 


<3000 


2000-3 xlO 5 


2xl0 4 


lxlO 4 


H 2 CO 


<3xl0 14 


<8 


4-19 


6.5 


3.9 


CH3OH 


<3xl0 14 


<7 


70-2000 2xl0" 4 


4xl0" 4 


SO 


<3xl0 13 


<0.7 


4-37 


0.05 


0.14 


HCN 


<4xl0 12 


<0.09 


2-18 


200 


500 


S0 2 


<3xl0 13 


<0.8 


4-19 3xl0" 4 


0.02 


CS 


<7xl0 12 


<0.2 


4-13 


0.17 


0.70 


CN 


<3xl0 13 


<0.7 


1.4-5 


2.6 


250 


HC 3 N 


<4xl0 12 


<0.1 


0.2-0.8 4xl0" 4 


9xl0" 4 


SiO 


<4xl0 12 


<0.09 


0.6-40 







Range of values got from the molecular outflows in LI 157 
(Bachiller & Perez Gutierrez \TW7l . NGC 1333 IRAS2A 
(J0rgensen et aL lSTOTl and BHR71 (Garay et al lTTOl 
Theoretical ratios derived for shocked molecular gas under 
a strong UV fields: see more details in § 15.41 



4.3. The High Velocity Region 

The column density for the HCO + and the upper limits 
(at 3a) for the rest of the molecules were estimated at the 
intensity peak of the high velocity HCO + emission (see 
Fig. El, which is ~ 3" northwest of the HH 2 knot E. 
The column densities are estimated from the integrated 
emission over the —5 to 5 and 11 to 21 km s" 1 velocity 
ranges, avoiding the contribution from the ambient gas of 
the Ahead Core (and the SO2 Clump), and assuming an 



excitation temperature of 30 K. Table HJ1 shows the column 
densities and the relative abundances of several molecules 
with respect to HCO + , and for comparison we also give 
typical relative abundances in shocked regions of molecu- 
lar outflows driven by low mass protostars. It is clear from 
this table that the HCO + abundance is strongly enhanced 
or, alternatively, that the other molecules are severely de- 
pleted in the High Velocity Region (see ? I5.4|> . 



12 



Girart et al.: Radiative and dynamical perturbations around HH 2 



□ Ahead Core 
1 S0 2 Clump: 14 

| S0 2 Clump: II 

West Core 



CO/10 3 H 2 CO HCO+ CH,OH SO HCN S0 2 CS CN 

Fig. 15. Comparison of the column densities of several 
molecular species for the Ahead Core, the West Core and 
the positions II and 14 of the SO2 clump. Position II is 
closer to HH 2 than 14 (see Fig. Upland Table El . Asterisks 
on the bar means upper limit of the column density. 



5. Discussion 

In this section, we attempt to provide possible explana- 
tions for the origin and structure of the observed molecu- 
lar emissions in the frame of UV illumination models. We 
do not in any way attempt to model any of the regions in 
great detail. From now on we will define the SO2 Clump 
as the molecular component within the Ahead Core traced 
by SO, SO2 and CH3OH and whose chemical properties 
were studied in Papers I and II (see Fig. 2 from Paper 
I). The SO2 Clump appears located in the region of the 
Ahead Core facing HH 2 (and thus is the molecular region 
in the Ahead Core more exposed to the UV radiation). 
Therefore, we split the discussion of the SO2 Clump and 
Ahead Core in two different subsections. Figure ITB1 shows 
the column densities for the Ahead Core, SO2 Clump and 
the West Core. 

5.1. The SO2 Clump 

In § !4.1.3l we found an SO excitation temperature and rel- 
ative molecular abundance enhancement within the SO2 
Clump towards HH 2 (see Figures IT21 and H3jl . The excita- 
tion temperature enhancement indicates a density increase 
towards HH 2, which may suggest compression from the 
wind. Although the position at which the molecular en- 
hancement occurs coincides very well with HH 2 knot E, 
the most likely source of UV radiation are the high exci- 
tation knots H and A (e. g., Bohm et al. 1993). This could 
be due to a projection effect: i. e., the molecular gas spa- 
tially coinciding with knot E could be the closest to knot 
H and A. 

It is possible that the photoelectric heating mecha- 
nism due to the strong UV radiation could cause an in- 
crease in temperature in the molecular gas. We investi- 



gated this scenario by running a PDR model that self- 
consistently computes the temperature at the edge of the 
clump (see Papadopoulos, Thi & Viti 2002 and Bell et al. 
2005) . We found that with a radiation field of ~40 Habing 
(Molinari & Noriega-Crespo HSH 1 Habing = 1.6xl0 -3 
erg cm~ 2 s^ 1 ) a temperature of 50 K is reached 1.0 mags 
into the clump for densities between 10 5 and 10 6 cm -3 . 

As already mentioned there is a definite gradient in 
both density and abundance. In light of this result, and 
recalling the conclusions of Paper II it is clear that the 
molecular enhancements found in the SO2 Clump do not 
come from a single density component gas. We have there- 
fore computed the observed column densities and the tem- 
peratures at four positions along the slice of Figure ^3 
see Table 0] In this table we also give the distance be- 
tween the positions selected and the HH 2 knot H (which 
together with knot A are the likely source of the UV ra- 
diation). We then ran multi-density components models, 
using the chemical model employed in Paper II and revised 
visual extinctions. The physical and chemical parameters 
of these models were as those in B2 and in B7 in Paper II, 
but we used a number density structure with a minimum 
and maximum number density of respectively 1x10 s and 
3xl0 5 cm~ 3 . The maximum number density is as in B2 
and B7 of Paper II, while the minimum density is deter- 
mined by applying the density law derived by Tafalla et 
al. ( 2002) for starless cores. Since in our models the visual 
extinction is a function of the column density and there- 
fore of the density, the visual extinction at any point in the 
modelled clump has changed from B2 and B7 accordingly. 

We find that it is possible to match, within a factor 
of 5, all the species everywhere if the visual extinction 
(along the line of sight to the source) is ~ 3 mags for all 
positions, and for ages < 1000 yrs (after it starts being 
irradiated). Since we know the distance of the emitting re- 
gion to the UV source (cLhh2 in Table 4) , a constant visual 
magnitude suggests one of the following three scenarios: 
(1) a decreasing (from HH 2) density structure with den- 
sities of 1.4xl0 5 cm~ 3 , 7.4xl0 4 cm" 3 , 5xl0 4 cm~ 3 and 
3.5xl0 4 cm -3 for positions II, 12, 13 and 14, respectively 
(see Table Q); (2) a clumpy structure, where the density 
and size of each clump is constant (and hence the column 
density within each clump along the line of the source), 
and the difference in the observed column densities is sim- 
ply due to the fact that as the radiation field crosses sev- 
eral clumps, it gets attenuated. This implies that the four 
positions correspond more or less to four clumps along the 
same line with respect to the source; (3) a third scenario 
may involve a more complicated picture where the line 
joining knots A, H and the SO2 Clump is not perpendic- 
ular to the line of sight, hence the radiation field may be 
impinging on the SO2 Clump from the top (or the bot- 
tom) along the line of sight to the observer. Probably the 
most likely scenario is a combination of the first and third 
scenarios. The first one agrees with the density increase 
towards HH 2 derived from the SO line ratios. The third 
one may agree with the distribution of the SO, SO2 and 
CH3OH with respect to the dust: these molecules appear 



Girart et al.: Radiative and dynamical perturbations around HH 2 



13 



located in the exposing face (with respect to HH 2) of the 
strongest 0.85 mm dust peak from the Ahead Core. In 
this case, the SO2 Clump may not be a true clump, but 
just the illuminated face of a slightly larger clump, whose 
center is located at the position of the 0.85 mm dust peak. 

On the basis of this model, we conclude that the HH 2 
is the source of 'illumination' and that the SO2 Clump 
must have substructure (a density increase towards HH 2, 
as a result of compression from the wind) within it (re- 
iterating a main conclusion of Paper II) and that this ac- 
counts for the chemical variation across the core. However, 
it is not possible to deduce from these studies whether the 
SO2 Clump is a core distinct from the Ahead Core or is 
simply a face of Ahead Core strongly affected by the ra- 
diation from HH2. 



5.2. The Ahead Core 

The Ahead Core clearly shows some molecular abundance 
differences with respect to a typical starless dense core: 
for example, comparing the Ahead Core abundances with 
those from L1544 and L1689B (J0rgensen, Schoier, & van 
DishoeckHU® i both CO and HCO+ are slightly enhanced 
(a factor ~ 2-8), whereas the CS is depleted by a factor 

9- 26. We now explore whether, as for the SO2 Clump, the 
abundances can be reproduced as a consequence of the 
UV field from the HH 2 shock impinging on the gas and 
dust of the core. 

Excluding the SO2 Clump, the Ahead Core has a large 
size, ~ 0.1 pc, hence it is reasonable to assume that the 
species we are detecting may all be tracing different den- 
sity components in this core. We used the models in Paper 
II to determine whether the chemistry of the core can be 
explained by a weak radiation field, because further away, 
that is impinging on a fairly dense core. We find that the 
results at early times (< 500 yrs) of models from both Grid 
A and B from Paper II may match the observations fairly 
well (within a factor of 5, and with a radiation field of ~ 

10- 20 Habing) if one assumes that i) different densities 
along the line of the source are present and that ii) dif- 
ferent molecules trace different components. For example, 
CS and CO emission seem to come mainly from the south- 
ern part of the Ahead Core, correlating well with a sec- 
ondary dust peak, while HCN, H2CO and HCO + trace the 
whole core. These considerations suggest that the Ahead 
Core may not be an homogeneous and roughly spherical 
core but it may be composed of small clumps similar to 
those found in L673 (Morata, Girart & Estalella 1271775) . 
These BIMA observations reveal small-scale clumps (typ- 
ically a few x 10~ 2 pc) in lines of HCO+, N 2 H+, and CS. 
However, the maps of these species did not coincide; this 
is probably a consequence of the time-dependence of the 
chemistry (Garrod et al. 2005). A detailed modelling of 
this region, taking into considerations these assumption 
should be attempted. 
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Fig. 16. Sketch of the West Core, the PDR-like structure, 
and shocked regions that lie behind the West Core with 
respect to the observer point of view. 

5.3. The West Core 

The well defined ring structure of optical and mid-IR 
emission (Warren-Smith & Scarott 1999, Lefloch et al. 
2005) around the West Core suggests that its geometry is 
quite spherical, with a radius of ~ 0.03 pc. The size of the 
West Core is similar to those from the transient clumps 
resolved in L673 (Morata et al. l2UTJ5|) . 

The presence of [SII] and Ha in the ring indicates that 
the West Core is already being shocked by the protostel- 
lar winds of the YSO VLA 1. But how can the wind from 
VLA 1 form a ring structure around the West Core in 
the plane of the sky? A possible answer is that the shock 
comes from the 'dark face' of the West Core, as shown 
in the sketch of Figure ITCl In this case, the cold and dy- 
namically unperturbed molecular gas and dust is prob- 
ably masking the shock where the column densities are 
high enough. The gas column density at the dust peak is 
N(H 2 ) ~lxl0 22 cm -2 (see § I4.2f> . which implies a vi- 
sual extinction of about 11 magnitudes. This is sufficient 
to create an apparent hole in the [SII] and Ha emission 
at the center of the West Core. Following standard ex- 
tinction laws (e.g., Mathis ll99U|) . the mid-IR extinction 
ranges between ~ 0.2 and 0.7 (i.e., an intensity decrease 
between 20% and 50%), which also can account for the 
observed mid-IR ring-like structures. The 13 CO, HCO + , 
HCN and CS spectra from the West Core show two lines, 
one at a i>lsr of ~ 6.3 km s _1 and quite narrow (the 
FWHM of the 13 CO 1-0 line is ~ 0.48 ± 0.03 km s" 1 ), 
and a second line, weaker but broader, Av ~ 1.5 km s _1 , 
centered at a vlsr of ~ 8.7 km s _1 (see Fig. 01 . The nar- 
row line possibly traces dynamically unperturbed molecu- 
lar gas from the West Core, whereas the broader line may 
trace a molecular layer of interaction between the West 
Core and the shocked material. Further observations are 
required to confirm this hypothesis. 

The ring-like mid-IR emission around the West Core is 
tracing warm gas interpreted by Lefloch et al. j2005| as a 
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PDR created by UV radiation field of 20-40 Habing. Thus, 
the properties of the West Core, i.e., a unperturbed dy- 
namically core surrounded by PDR, suggest that the West 
Core is in a similar situation as the SO2 Clump. However, 
the chemical composition of the West Core seems differ- 
ent: the column densities for most molecules are a factor 
of ~ 10 lower in the West Core (Table EJ) than in position 
14 (Table©, except for HCN, CS and CN which are of the 
same order. 

The major constraints for the unperturbed cold molec- 
ular gas from a chemical point of view are that i) the 
line widths are narrow and hence the gas is probably not 
shocked; ii) there is an enhancement of HCO + with re- 
spect to a typical dense core but no other enhancement 
is observed and in fact CH3OH and SO are not even de- 
tected. The BIMA observations indicate a quiescent, pos- 
sibly cold (< 30K) clump. On the other hand, Lefioch et 
al. PIH^I detect warm gas interpreted as a PDR at the 
edge of this clump. We have run a PDR model (with a 
radiation field ranging from 20 to 40 Habing: Lefioch et 
al. 2005) of this region to estimate whether the presence 
of a PDR at the edge of the West Core can coexist with 
the quiescent dense gas we observe. We found that as long 
as the average density of the West Core is between 5xl0 4 
and 5xl0 5 cm -3 , then it is possible for the central re- 
gion of the core to be cold. This range of densities are in 
agreement with those derived from the PDR (Lefioch et 
al. 12005(1 . In particular we find that for a density of 10 5 
cm~ 3 , the gas temperature reaches 42 K at 1 mag and 27 
K at 3 mags. 

By combining the abundances of the modeled PDR 
with those from the models in Paper II we find that the 
best match is reached by a 2-component structure made of 
a fairly old quiescent, photo-processed clump (^1000 yrs 
of age) of thickness less than ~ 2 mags plus a PDR com- 
ponent; such a combination gives the right abundances of 
CO, CS, H 2 CO and CH 3 OH (coming from the old dense 
clump), and of HCO + (coming mainly from the PDR); we 
find that the PDR needs to be at least 1 mags. 

An alternative view is that the mid-IR emission comes 
not from a conventional PDR but from a warm interface 
with a PDR created in the interaction of the winds of 
the YSO VLA1 with the West Core. Nguyen, Hartquist & 
Williams ( 20(3} have explored the physics and chemistry 
occurring in these interfaces, and have shown that the 
thermal pressures in the interfaces should significantly ex- 
ceed those within the cloud. Resulting temperatures are 
expected to be ~ 10 3 K and dust emissions may therefore 
mimic those of a PDR. The results of Nguyen et al. ( 2001 ) 
suggest that, if the interface occupies a few percent of 
the core, then some species may achieve detectable abun- 
dances. While CO and CS are not particularly enhanced 
above cold cloud abundances, in some conditions HCO + 
may be enhanced by up to two orders of magnitude above 
its cold cloud abundance. H 2 CO appears to be somewhat 
enhanced; results for methanol are not reported. Thus, 
it appears possible that the dynamical interaction of the 
wind and the clump may lead to a characteristic chem- 



istry in the interface. Such a model requires more detailed 
study, perhaps along the lines of that by Lim, Rawlings & 
Williams ifTMTjl . 

5.4. The High Velocity Region 

The properties of the high velocity HCO + emission as- 
sociated with HH 2 are similar to those found in the 
NGC 2071 outflow (Girart et al. ITMUjl : there is a good 
correlation with the shock-excited H 2 (see Fig. |SJ) and 
the emission has a monotonic enhancement with the flow 
velocity (Dent et al. 127703) . The high velocity HCO+ emis- 
sion is also well correlated with optical knots with small 
proper motions (e. g., L, E), with (tangential) velocities in 
the 20-80 km s" 1 range (Bally et al. 1271172) . Since HCO+ 
is expected to be strongly enhanced in turbulent mixing 
layers associated with low-velocity (40 km s _1 ) shocks 
(Taylor & R aga IMTH|> . Dent et al. pT(T3) suggested that 
this flow velocity enhancement dependence could be ex- 
plained if the HCO + is formed at high shock velocities 
and then gradually mixes with ambient unenhanced ma- 
terial. This is like the wake chemistry explored by Lim et 
al. (1999) and implies that a UV radiation may not be 
necessary in the chemical processing of the HCO + (Dent 
et al. 12003}! . However, in the mixing layer scenario one 
would expect other molecules, such as CS, H 2 CO, SO, to 
be also enhanced (e. g., Viti, Natarayan & Williams 2002), 
but these are undetected in the High Velocity Region. 

Analysis from the mid-IR H 2 line emission shows that 
the 0-0 S(n) lines have contributions from hot (~ 1000 K) 
and warm (~ 300 K) components (Lefioch et al. 2003). 
The hot component arises from compact regions (Ss 2"), 
with shock velocities of 20-30 km s . The warm com- 
ponent arises from more extended and slower shocks (10- 
15 km s _1 ) in a denser medium, ~ 10 5 cm~ 3 . The fact 
that the high velocity HCO + emission is better correlated 
with the S(2) H 2 line than with higher excitation H 2 lines 
and that this line traces the warm H 2 component sug- 
gest that the high velocity HCO + probably arises from or 
around the H 2 warm component. 

As for the other three regions, we have attempted to 
qualitatively model the high velocity gas. The chemical 
anomaly of this gas is that it differs substantially from a 
typical shocked region (see Table 0). In addition, we have 
tried to see whether UV radiation alone can reproduce 
the HCO + overabundance and found that this is not the 
case. We have investigated whether it is possible, via a 
combination of heating (e.g. induced by the passage of 
a low velocity, non-dissociative shock, as the HCO + and 
warm H 2 correlation suggests) and UV chemistry, for a gas 
to have such a composition as we find here, i.e. abundant 
HCO + without the presence of other shock tracers such 
as SiO and CH 3 OH. 

The models used for this region are similar to those 
employed in Viti et al. ( 2004 ) where we explore the chem- 
ical evolution of low velocity, chemically rich clumps ob- 
served along the main axis of chemically rich outflows. The 
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chemical model simulates the clump formation (by free- 
fall collapse) and its subsequent interaction with the out- 
flow; here we only consider the scenario where the density 
structure formed before the advent of the outflow (see Viti 
et al. 120041 for a more detailed explanation of the model) . 
In addition, here we also include the presence of a strong 
radiation field. After a first stage where densities of up to 
10 5 cm -3 are reached, the region is heated and affected 
by a very strong UV field. We simulate a short period 
(~ 100 yrs) of high temperature (~ 1000 K), followed by 
cooling, down to a temperature of 200 K. We employ a UV 
field of 100 Habing: this value is derived from the consider- 
ation that if the edge of the West Core is affected by a UV 
field of 20-40 Habing, then since the High Velocity Region 
is closer to the HH shock by a factor two the radiation 
field impinging on it should be about four times stronger. 
TableElreports the column densities computed for species 
observed in this shocked gas. Observations show that the 
only detected molecule is HCO + while for the other species 
we have upper limits. We tested two models: Model A 
where we assumed a core size of ~ 0.02 pc that, at a den- 
sity of 10 5 cm -3 , corresponds to a visual magnitude of 
~ 4.8 mags; and Model B, as Model A but for a smaller 
clump of 0.01 pc, corresponding to a visual magnitude of 
2.8 mags. We ran the chemistry for 100,000 yrs. 

For model A we found that the best match with ob- 
servations is reached at 40,000 yrs, where HCO + peaks 
with a fractional abundance of 5 x 10~ 9 ; the gas has by 
then completely cooled down (to 200 K). For model B we 
found that HCO+ peaks earlier, at 1000 yrs; the HCO+ 
fractional abundance is 6 x 10 -9 . For both models, we 
report in TableElthe ratios for selected species to be com- 
pared with those observed and shown in the same table. 
If we consider that the HH jet is traveling at velocities of 
hundreds of km s -1 , then model B seems more realistic 
because of the shorter timescales involved. 

There is a remarkable correspondence between the 
computed peak abundance of HCO + , and the observed 
column density: the latter is derived from observations 
to be 4.1 x 10 13 cm -2 which implies that for a 0.01 pc 
clump (Models B) the observed fractional abundance is 
- 8.9 x 10" 9 . 

In model B, HCO + reaches ~ 10 -9 very early, at 
about 500 years; it then increases but after about 1500 yrs 
starts being destroyed (a consequence of the lower Ay). 
However, during these 1000 yrs, which include the cooled 
phase, HCO + is ~ 5xl0 -9 for about 400 yrs and dur- 
ing that period the ratios of most of the other species to 
HCO + are consistent with the observed ratios (CO/HCO + 
is higher than the ratio estimated from the observations by 
Dent et al.[2003). In Table we report theoretical ratios 
for all the species, (but not for SiO, which is not included 
in our models), for which either a detection or an upper 
limit is derived from observations. We chose to show the 
ratios at the time when HCO + peaks in our models. The 
species that differ most from observations are HCN and 
CN. As extensively explained in Paper II, there is a funda- 



mental problem with these two species that at the moment 
we are unable to explain. 

In summary, we have found two suitable scenarios that 
match the observations of the High Velocity Region except 
for the CO, CN, and HCN abundances. The essential fea- 
ture is of a warm chemistry (possibly due to the passage of 
a low velocity, non-dissociative shock) which is subjected 
to stronger than ambient UV field, a scenario not previ- 
ously explored. 

5.5. Is the molecular gas ahead of HH 2 being driven 
out by the VLA 1 winds? 

The kinematical properties described in § 13 .41 suggest that 
the Ahead Core is being driven by the protostellar winds 
from VLA 1 . If the velocity gradient observed is due to the 
swept-up effect from the winds, then these winds should 
have a low collimation since the dense molecular struc- 
ture traced by the HCO + (including the Ahead Core) is 
quite extended. Low collimated winds sweeping up dense 
molecular gas have being observed around the HH 34 sys- 
tem (Anglada et al. 1155511 and XZ Tauri (Welch et al. 

nam . 

We considered a model of a spherically symmetric stel- 
lar wind, with constant velocity V w and mass-loss rate 
M w , sweeping up the ambient material as a snow-plow 
and accumulating it into a shell. Anglada et al. (1995) 
obtains analytical solutions in the case of non-negligible 
ambient cloud pressure pAV^, for ambient cloud density 
power-law distributions, p ex r~ Q , with a — 2 and a — 
(constant density) . A V c is the mean squared velocity (tur- 
bulent plus thermal) in the ambient cloud. The parame- 
ters of the shell used were those observed for the Ahead 
Core, i?shcii = 5xl0 4 AU, the distance from the Ahead 
Core to VLA 1, 14hcii = 3.0 km s~ x , the maximum ve- 
locity relative to the systemic velocity observed in the 
Ahead Core, and A/ s i lc ii = 3.8 Mq, the mass derived from 
the dust (Dent et al. [2003 ), and a mean-squared veloc- 
ity AV C = 0.6 km s _1 , the typical observed line width. 
Table shows the physical parameters obtained for the 
density of the ambient cloud, the age of the shell, the mo- 
mentum rate of the stellar wind, and the mass-loss rate 
for a wind velocity of 100 km s _1 . 

In our case we are only observing the southern part 
of the region surrounding the powering source of the HH 
1-2 outflow, VLA 1, and the Ahead Core subtends only 
a small solid angle fl < 4ir sr around VLA 1. However, 
the results obtained for a spherically symmetric wind are 
applicable to the case of a conical stellar wind. In general, 
if the solid angle of the wind cone Q w is greater than the 
shell solid angle fi s heii> the momentum rate, and mass-loss 
rate obtained in the spherically symmetric case have to be 
increased by a factor Oui/^shell- &w is at least 4 times the 
observed fishe!!! a factor of 2 comes from the fact that 
the wind is bipolar, and an additional factor of 2 comes 
from the fact that the outflow axis is very close to the 
plane of the sky (Solf & Bohm 1991) and that the Ahead 
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Table 7. Physics of the wind driven Ahead Core a 



Parameter Value 








2 


tsheii (yr) 


3.9xl0 4 


7.9 xlO 4 


n(H a ) ( cm~ 3 ) 


1.3xl0 3 


4.3 xlO 2 c 


P w (Mq yi^ 1 km s" 1 ) 
Mt (Mq yr- 1 ) 


2.9xl0 -4 


1.5xl0~ 4 


2.9xl0 -6 


1.5xl0~ 6 



Obtained for shell parameters -Rsheii = 5xl0 4 AU, Vghell = 

3.0 km s" 1 , Msheii = 3.8 Mq, and AV C = 0.6 km s" 1 . 
b Power index of the the radial density distribution, n(Eb) oc 

r~ a 

c Value for a distance of 5xl0 4 AU. 

d Value obtained assuming a wind velocity of 100 km s _1 . 



Core lies in the foreground face of the outflow lobe (see 
? I3.4|I . and, for symmetry, the wind has to exist also in the 
background face. In conclusion, the density, momentum 
rate, and mass-loss rate obtained in Table should be at 
least a factor of ~ 4 higher. 

As the outflow axis lies near the plane of the sky, pro- 
jection effects have to be accounted for. The de-projected 
velocity and distance are 14heii/ sini and -R s heii/ cosi, re- 
spectively, where i is the inclination with respect to the 
plane of the sky. However, the inclination is the biggest 
uncertainty. The dynamical age of the HH 1-2 outflow, 
~ lxlO 4 yr ago fQgura H995f) and the momentum rate of 
the HH 1-2 jet is 8xl0" 4 M & yr" 1 km s" 1 (Chernin & 
Masson 1995} ■ Given the uncertainties, the wind driven 
model can roughly account for these values for inclina- 
tions in the range 30°-45°. For such a range t s heii does 
not change much (a factor 0.6-1) and the outflow mo- 
mentum rate (for a = 2) becomes ~ 1-2 xlO -3 Mq yr -1 
km s . 

The total mass of the molecular gas and dust in the 
HH 1-2 region, derived from ammonia observations, is ~ 
52 Mq (Torrelles et al. ©94) . Taking into account the 
mass of Ahead Core and of the West Core (~ 4.1 Mq), 
the fraction of the molecular mass that is being driven 
out by the VLA 1 winds in the HH 2 lobe is ~ 8 %. Since 
it is possible that a similar situation may happen in the 
northern lobe around HH 1 (Torrelles et al. ©93 found 
also 'quiescent' ammonia clumps ahead of HH 1) the total 
mass driven by the powerful VLA 1 is likely 10-15 % of 
the total mass in the HH 1-2 region. 

Finally, it is interesting that while the molecular abun- 
dances may be explained by chemical models of pre- 
existing quiescent and transient clumps irradiated by a 
strong UV field (see previous sections and Paper II), the 
kinematics of the regions suggests that the molecular gas 
is indeed driven out by protostellar winds and the SO line 
ratios in the SO2 clump indicate a density increase towards 
HH 2. These results seem to be apparently in contradic- 
tion. Although beyond the scope of this paper and possi- 



bly quite speculative, it is possible that these clumps have 
formed due compression by winds but in a somewhat sim- 
ilar way (though probably faster) to the formation of the 
transient clumps. The presence of transient clumps within 
molecular clouds has been reported recently both observa- 
tionally (Morata et al. l2()()5|l and theoretically (e. g., Falle 
& Hartauist 120021 Vazquez-Semadeni et al. 2005). 

6. Summary and Conclusions 

We have carried out an extensive observational study 
(from BIMA data) and attempted a preliminary theoret- 
ical investigation of the molecular gas around HH 2. The 
BIMA maps show a very complex morphological, kine- 
matical and chemical structure of the molecular gas. For 
clarity we divided the observed region in four subregions 
taking into account the properties of the different molec- 
ular species and the properties observed at other wave- 
lengths: 

— The Ahead Core, located ahead of HH 2, which has a 
size of ~ 0.1 pc and a mass of 3.8 Mq (Dent et al. 
2003). The observed molecular abundances from the 
Ahead Core differ from the typical values of low-mass 
protostellar envelopes. A weak UV field (weak because 
attenuated when passing through the SO2 clump) orig- 
inated in HH 2 can account for the observed values if 
different species trace different layers. The differences 
between species within the Ahead Core suggest that it 
is not homogeneous but probably composed by small 
clumps. 

— The SO2 Clump is a molecular component within the 
Ahead Core that is more exposed to the UV radia- 
tion from HH 2. The analysis of the observed molec- 
ular lines indicates an increase of density and relative 
molecular abundances towards HH 2. The photoelec- 
tric heating produced by the strong UV field could 
probably create a temperature enhancement at the 
edge of the clump. A four point chemical analysis of 
the SO2 Clump confirms the conclusion of Paper II 
that the clump must have substructure within it if its 
chemistry is due to the HH 2 UV illumination. In par- 
ticular, an increase of density towards HH 2 may ac- 
count for the chemical properties, which would suggest 
compression from the VLA 1 winds. From the overall 
properties of the SO2 Clump we cannot say whether it 
is a clump distinct from the Ahead Core, or, whether it 
is part of the face of the Ahead Core exposed to HH 2. 

— The West Core is a molecular structure with a radius 
of 0.03 pc and a mass of ~ 0.3 Mq surrounded by 
a ring structure of shocked gas traced by [SII] and 
Ha (Warren-Smith & Scarott ITMty and of a PDR- 
like structure traced by mid-IR hot dust and PAH 
emission (Lefloch et al. 120050 . This ring-like structure 
is likely not real but a consequence of the fact that 
the West Core is in the foreground with respect to 
the shocked and hot component. The properties of the 
West Core and the SO2 Clump are somewhat similar 
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(i. e., apparently quiescent gas irradiated by a strong 
UV field), although their chemistry is different. We find 
that the chemistry of the West Core can be best ex- 
plained as arising from a combination of an old photo- 
processed dense clump, responsible for the emission of 
CO, CS, H 2 CO and CH3OH, and a PDR, from where 
most HCO + emission comes. Alternatively, the mid- 
IR emission may not come from a conventional PDR 
but from a warm interface with PDR created in the 
interaction of the VLA 1 outflow with the West Core, 
which may also account for the HCO + emission. 
— The High Velocity Region, associated with HH 2, is 
traced by the HCO + but not by other molecular shock 
tracers such as SiO, CS and CH3OH. The presence 
of enhanced HCO + and the lack of other shock trac- 
ers can be accounted by the interaction of the VLA 1 
outflow with a pre-existing dense clump via heating 
(possibly due to the passage of a low velocity, non- 
dissociative shock) and by the presence of strong UV 
field (coming from the nearby high excitation HH 2 
knots), although our models predict excessive abun- 
dances of HCN and CN. 

The overall main conclusion of this work confirms the 
findings of Paper I and II, by demonstrating that in ad- 
dition to the strong photochemical effects caused by pen- 
etration of the UV photons from HH 2 into molecular 
cloud, a range of complex radiative and dynamical inter- 
actions occur. These generate characteristic chemical sig- 
natures for each type of interaction. Thus, despite the ap- 
parent 'quiescent' nature of the molecular cloud ahead of 
HH 2, the kinematical properties observed within the field 
of view suggest that the cloud is not 'quiescent' but it is 
possibly being driven out by the powerful winds from the 
VLA 1 protostar. If so, this would imply that a significant 
fraction (10-15 %) of the total mass in the HH 1-2 region 
is being disrupted from the original molecular dense core 
around the protostar. Given the large size of the molecular 
emission, these winds should be of lower collimation that 
the jet associated and responsible of HH 1 and 2. 
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Fig. 1. Channel velocity (for Av = 1 km s" 1 ) BIMA contour maps of HCO+, H 13 CO+, H 2 CO, HCN and DCO+ 
(the maps are not corrected for primary beam attenuation). For the HCO + the contours are -4, -3, 3, 4, 6, 9, 12, 15, 
18, 21, 24, 27, 30, . . .80 times the rms noise of the map, 62 mJy beam" 1 . Contours are -4, -3, 3, 4, 5, 7, 9, 11, 
31 times the rms noise of the maps: 55 (H 13 CO+), 120 (H 2 CO), 36 (HCN), and 140 mJy beam" 1 (DCO+). The CN 
iV=l-0 map was obtained by combining the following hyperfme transitions: F=5/2-3/2, J=3/2— 1/2 F=3/2-l/2 and 
,7 = 1/2-1/2 F=3/2-3/2 (note that Table [21 give the rms of the maps previous to this combination procedure). The 
synthesized beam is shown in the lower left corner of the right panels. The grey ellipsoid shows the position of the 
brightest HH 2 knots (A, B, D, H). The two filled triangles show the position of HH 2 knots E (northern one) and L 
(southern one). The filled square shows the position of the driving source of HH 2. 
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Fig. 2. Channel velocity BIMA contour maps of three SO lines, SO2 and CH3OH (the maps are not corrected for 
primary beam attenuation). Contours are -4, -3, 3, 4, 5, 7, 9, 11, ... 25 times the rms noise of the maps: 95 (SO 
3 2 -2i), 85 (SO 22-lx), 73 (SO 2 3 -l 2 ), 33 (S0 2 ), and 50 mJy beam^ 1 (CH 3 OH). The CH3OH map was obtained 
by combining the 2o-loA + and 2_i-l_i E transitions (note that Table El give the rms of the maps previous to this 
combination procedure). 




Fig. 3. Channel velocity BIMA contour maps of 13 CO, C 18 0, CS, CN and C3H2 (the maps are not corrected for 
primary beam attenuation, except for 13 C0). Contours are -4, -3, 3, 4, 5, 7, 9, 11, ... 25 times the rms noise of the 
maps: 85 ( 13 CO), 75 (C 18 0), 83 (CS), 50 (CN), and 52 mJy beam" 1 (C 3 H 2 ). The HCN J=l-0 map was obtained by 
combining the three hyperfme transitions: F=2-l, F=l-1 and F=0-1 (note that Table give the rms of the maps 
previous to this combination procedure). The C 18 1-0 maps were obtained combing 109.7-113.0 and 109.9-113.4 GHz 
(LSB-USB) frequency setups. 



